T 




00 

o 



< 

o 


•—J 

CD 

cn 

o 



lr- .... - . 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



TECHNICAL MEMORANDUM 1408 


THE PRINCIPLES OF TURBULENT HEAT TRANSFER 

By H. Reichardt 


i 

! 

* 


Translation 

“Die Grundlagen des Turbulenten Warmeliberganges. ,, Archiv £. die 
gesamte Warmetechnik, no. 6/7, 1951. 





*• 



Washington 


September 1957 


AFMDC 

TECHNICAL LIBRARY 
AFL 2811 


TECH LIBRARY KAFB, NM 


»CK-1 


MCA TM 1408 



MTIOML ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1408 
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1. INTRODUCTION 

The literature on turbulent heat transfer has in the course of years 
attained a considerable volume. Since this very complicated problem has 
not as yet found a complete solution, further studies in this field may 
be expected. The heat engineer must therefore accommodate himself to a 
constantly increasing number of theories and formulas. Since the theo- 
ries generally start from hypothetical assumptions, and since they con- 
tain true and false assertions, verified knowledge and pure suppositions 
often being intermingled in a manner difficult to tell them apart, the 
specialist has difficulty in forming a correct evaluation of the indi- 
vidual studies . '--•■■■■ 

The need therefore arises for a presentation of the problem of tur- 
bulent heat transfer which is not initially bound by hypothetical as- 
sumptions and in which the already known and that which is still unin- 
vestigated can be clearly distinguished from each other. Such a 
presentation will be given in the present treatment. 

The following brief remarks may be made with regard to the develop- 
ment of the theory of local heat transfer. 1 The first to recognize 
the int ima te relation between heat transfer and flow resistance was 
0. Reynolds (ref. 17). The considerations of Reynolds hold, however, 
only for fluids with special properties (according to present termin- 
ology, they are fluids whose Prandtl number is l). 


* w Die Grundlagen des Turbulenten Warmeuberganges . r T Archiv f . die 
gesamte Warmetechnik, no. 6/7, 1951, pp. 129-142. 

^This is taken to mean the direction of investigation which has for 
its object the derivation of the law of heat transfer from local flow 
► processes . In contrast to this is the semiempirical method of similar- 

ity considerations as developed by W. Nusselt. This model investigation 
dispenses with the knowledge of the individual processes. However, for 
this reason it is in the position to supply simple practical formulas 
even in complicated cases . 
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The first practical useful formula for turbulent heat transfer was 
derived in 1910 by L. Prandtl (ref. 12). Prandtl, like Reynolds, started' 
from the assumption that heat and momentum are transferred by the same v 

turbulence mechanism. For simplifying the computation, the friction 
layer was divided into two sections : the turbulent boundary layer itself, 

and a thin layer close to the wall, whose flow was assumed as completely 
laminar . 

The Prandtl formula was found quite reliable for the representation 
of heat transfer of fluids with small Prandtl numbers. On the other 
hand, it is impossible to represent the measured heat-transfer coeffi- 
cients of fluids with high Prandtl numbers by the Prandtl formula (a 
correction is assumed which is adjusted to the experimental data) . 

The reason for this discrepancy for high Prandtl numbers lies in 
the too greatly idealized flow relations in the friction layer. Actually 
there is, of course, no completely turbulent boundary layer and also no 
completely laminar wall layer, but a continuous transition from the tur- 
bulent flow to the viscous flow in the Immediate neighborhood of the 
wall. 

By dividing the friction layer into three sections (namely, a tur- 
bulent region, a laminar- turbulent transition zone, and a laminar wall 
layer), as was done by Th. v. Karma'n (ref. 7) and H. Reichardt (ref. 14), 
a considerable improvement in the theory could therefore be attained . 

There still remained, however, certain contradictions which resulted 
from the conception of a completely laminar layer at the. wall . 

In order to remove these discrepancies, the author, in a later pa- 
per (ref. 15) , assumed an entirely continuous decrease of the turbulent 
exchange to zero near the wall. The results of this computation are 
contained in the present paper. In addition, still more general formu- 
las are given that form the basis for a later investigation of the effect 
of the temperature dependence of the constants of the material on the 
heat transfer. Also considered in this paper is the turbulent heat 
transfer for extremely ‘low Prandtl numbers. 2 


2. DERIVATION OF GENERAL FORMULAS 
(a) Definitions and Assumed Expressions 

The considerations will be restricted to those cases for which no 
unsteady changes of temperature occur at the wall to which the heat (or 
cold) is transferred from the flow. The wall is furthermore assumed as 
smooth . 


2 At the suggestion of Mr. B. Koch. 
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If the temperature of the wall is constant or its spatial changes 
are only slight, the heat flow and momentum flow at each point x have 
practically the same direction as y (or -y) perpendicular to the wall. 

The transfer of the momentum and heat is, on the one hand, effected 
through the molecular motion, and on the other hand, through turbulent 
exchange. The parts determined by molecular transfer will be denoted by 
the subscript m and those produced by turbulence denoted by the 
subscript t. For the flow densities of the momentum T and the heat 
q we therefore write 

T ■ T m + T t = T m (l + T t /c m ) (l) 


<1 - + It - + ‘ItAm) 


( 2 ) 


The following equations apply for the individual partial flows : 


and 


T^ m - h du/dy 

(3) 

r 0 - p 0 (du/dy) Q 

(3a) 

<3m » X dT/dy 

(4) 

qo = X Q (dT/dy) o 

(4a) 

r t = Ap du/dy 

(5) 

It - G p A q dT/dy 

(6) 


Through these formulas the coefficient of viscosity p and of heat 
conduction X, as well as the exchange magnitudes^ for the momentum 

and for the heat Aq, are defined (u and T denoting, respectively, 
the time mean value of ‘the flow velocity and the temperature at the dis- 
tance y from the wall; Cp is the specific heat). 

The subscript 0 refers to the wall. Since in the Immediate neigh- 
borhood of the wall only molecular transfers are possible, *to is identi- 
cal with 0 %)q and qo is identical with (<3m) 0 > while (t-^q and (qt) Q 

do not exist. 


^The exchange magnitude A T can be denoted as the turbulent friction 
coefficient, the magnitude CpAq as the turbulent heat -transfer 
coefficient. 
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The following equations are obtained from equations (5) and (3) 
and from equations (6) and (4), respectively: 

^tAm = A (?) 

a t Am - Pr* T tAm ■ P** AtA ( q ) 

where Pr* = Pr Aq/A T denotes the "generalized Prandtl number" (ref. 
14). Further , from equations (2) and (8) there is obtained 


Om = 1 

q 1 + Pr* Ar/p 


( 9 ) 


In order to be able to represent the temperature distribution and 
the heat transfer as functions of the flow magnitudes, still another re- 
lation is required that connects the dimensionless heat flow qAo w ith 

the dimensionless momentum flow t/t_. Since these two flows do not 
deviate strongly from each other in the neighborhood of the wall (where 
the main part of the heat transfer occurs) - at the wall itself we ac- 
tually have t /tq = aAo ~ -*• “ the following .expression, suggests itself: 

q/qo - (1 + k)T / t 0 _ (10) 

This equation defines a magnitude k which is small in the neighborhood 
of the wall and may therefore, to a first approximation/be neglected. 

The expression (10) was found to give reliable values for computa- 
tions involving medium and high Prandtl numbers (refs. 14 : and 15). It 
can, however, as will appear below, be successfully employed also for 
6mall Prandtl numbers (for which k is no longer small as compared with 
1 ). 


From the above definitions and assumed expressions there will now 
be derived the required equations for the temperature distribution and 
the heat transfer. The flow magnitudes A T /p, Pr*, and t'/tq are here 
to be considered as "given magnitudes." In order to assure the general 
validity of the formulas to be derived , no special assumptions in regard 
to the flow magnitudes will be introduced in this section. 


(b) Formulas for Temperature Distribution 
From equations (2) , (4) , (4a) , and (8) there follows 

to KDo ■ X ■! {1 + *** 


(11) 
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It is of advantage to introduce a dimensionless temperature difference 



which, is referred to the maximum temperature difference 0 at the dis- 
tance y = r (hence, 0 < 6 < 1 for 0 <.y/r< l) . For the temperature 

gradient r j , there is obtained from equation (ll) 


d-8 ( d6 \ ^0 

d(y/r) \<i(y/r)J 0 X 1 + Pi* Ar/u 


If equation (10) is used, there is further obtained from the above 



*•0 (1 + kH/Tp 
T 1 + Pr + A^/n 


a (y/r) 


(13) 


This equation represents the dimensionless temperature as a function 
of the dimensionless distance from the wall y/r. The factor of propor- 
tionality ^ is obtained from the condition that for y/r = 1 • 

we have •6 = 1. 

A first approximation of 6 is obtained by setting k equal to 
zero, X equal to Xq, and A^ equal to A^-. Under the integral there 

then appear only the magnitudes t/t o, A^/u, and Pr, considered as ini- 
tially given. 

For computing a second approximation of 6 , an equation is required 
for the heat flow appearing in the numerator of equation (13) . This 
equation is obtained from the continuity condition for the heat in the 
case under consideration. For example, in the case of the fully devel- 
oped flows in a pipe or in a plane rectangular channel, the following 
relations hold: 

py/r /P 1 

* cp ( 1 " r) d(y / r yj ^ cp ( 1 " d ( y / r ) 


i ( 1 - A = 1 - 
<10 V r / 


o 


(pipe) (14) 



I 


6 


MCA TM 1408 


-1 = 1_ 
10 


r>y/r /r 1 

d<P d(y/r )y 

0 /u 


d<p d(y/r) (channel) (15) 


where (p = u/U (ref. 14, p. 316). 


Substituted here will be the value of d of the first approximation 
of equation (13) . The approximation of q/q q, thereby obtained, yields 
a sufficiently accurate computation of the temperature distribution and 
of the heat transfer. 


Formula (13) is particularly adapted for computations with small 
Prandtl numbers . For high Prandtl numbers for which the integrand be- 
comes very small at even small distances y/r , it is more convenient to 
represent the temperature (at least as a first step) as a function of 
the velocity. 

For such representation there is required the velocity as a function 
of the distance from the wall. This relation is to be obtained from 
equations (l) , (3), (3a), and (7): 

- n (1 + (16) 

t o Wo 


After introducing a dimensionless velocity <p = u/U, which is referred 
to the maximum velocity at the edge of the boundary layer (at the dis- 
tance y/r), the following equation is obtained for the velocity gradient: 


d<p m f d(p > VO t Aq 

^ (y/r) Wfi ) 0 H 1 


(17) 


Combining equations (12) and (17) and employing equations (io) yield 


d 


/ddA ** C P0 (1 +a cM( 1 +k) 

U <pj 0 Pr 0 c 1 + Er* Vm 
C/Q 


(18) 


The proportionality factor (dd/d<p) 0 is obtained from the condition that 

for <p = 1 we likewise have d => 1. As regards the practical carrying 
out of the computations, the procedure with equation (18) is exactly the 
same as with equation (13) (see above). 


For the particular case Pr* = 1, equation (18) to a first approxi- 
mation gives the equality of d and <p. 


*r 


* 


* 




u 


T/PC 
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(c) Formulas for Heat Transfer 

It has become customary to express the heat transfer by dimension- 
less parameters. Generally used is the so-called Nusselt number., which 
is defined by the formula 


Hu = 


q0<3 

*0® U 


(19) 


where d denotes a characteristic length, (for example, the diameter of 
a pipe) and AT U is the temperature difference between the wall and a 
so-called mean flow temperature , 


of 


AT u = f ^ 


ATudf 


which refers to the quantity of heat f flowing through the flow cross 
section ^u m = ^ udfj . 

Since in the present paper all temperature differences are referred 
to the maximum temperature difference 0 at the distance y = r, it is 
convenient to introduce a dimensionless temperature difference, 

*u ■ AT u/® 

Setting further d = 2r and qQ = X^CdT/d y)o, there follows from equa- 
tion (19) 



( 20 ) 



formula 


follows from equation (13) for & = y/r 
is then obtained for the Husselt number: 


1. The following 


Nu« u 



*-o (1 + k HAo 

T 1 + Pr*A 


d(y/r) 


( 21 ) 


This formula holds for an arbitrary boundary- layer flow- In the 
case of laminar flow, A T - 0, so that the denominator of the integrand 

is equal to 1. For turbulent flow also, this denominator can at 
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least be set approximately equal to 1 if the Prandtl number is extremely 
low and A^/n is not too large (small Re number) . The quantities q/qQ 

and k are naturally different for the laminar flow than for the turbu- 
lent flow for vanishing Pr*A x /fi (fig. 5) . 

The applicability of formula (21) is problematical, however, since 
the variation of the heat flow q/qQ, for the case under consideration, 
must be known; here again a knowledge of the temperature distribution 
is required. Practically, therefore, the determination of Nu4 u is 
based on the somewhat inconvenient computation of the temperature dis- 
tribution which was presented in a formal manner in the preceding 
section. 

In order to be able to compute the heat transfer at the wall di- 
rectly, a formula is therefore required in which only known or easily 
determinable magnitudes appear. Such a formula could be obtained through 
a transformation of equation (2l) if it were possible' to introduce the 
factor Pr* into the numerator of the integrand. This, however, is un- 
fortunately impossible, because Pr* must be considered as a function 
of the distance from the wall through the magnitude A contained in 

it . Still another formula for Hu will therefore be derived in which 
the desired factor Pr* appears in the numerator of the integrand. 


For this purpose, equation (12) is multiplied by 


i2- ts, A i 


-P0 


a 


substituting equation (10) and integrating, there is then obtained 


p»4 


-PqA ' 1 


d4 


( <3* , 




ny/ r 


o 


U 0 Pr* (1 + k^Ao 
IjPro 1 + Pr*Ar/(a 


( 22 ) 


For 4 = y/r = 1 this is a determining equation for the temperature gra- 
dient to be substituted in equation (20) . Furthermore, using 


the notation^ 


(V) 

Ut 4 


^i 

C P^ L d4 

o Cp o AT 


(23) 


^Considered here is the mean value ^ 
temperature range 0 < 4 < 1. 


c Pq At 


d4 for the total- 
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yields 


Ha« u 



Pr»(l + k)T/x 0 mo d(y/r) 


1 + 


Pr^'/tx 


(21a) 


This equation, in spite of its different form, agrees in content with 
equation (21); if Ag/A T is considered as a constant, equation (21a) 

transforms directly into equation (21). 

The required transformation can now be made in equation (21a) by 
making the following substitution: 


Pr* 1 A Pc* - 1 

1 + Pr*A T /n 1 + A t /n (1 + Pr^/n) (1 + A T /ji) 

Making use of equation (17) and setting^ 

ol 


t = 


Pr*(l + A^/v) 

1 + Pr*A r /u kd<P 


k&p 




yields , from equation (21a) , 

dtp 


H U, U = *«(, 



(Pr* - l)d(p . „ 
1 + ?r*A T /p 


(24a) 


(25) 


(26) 


The form of this equation may further be improved by collecting 
factors together. Since 


2rPr 0 (dcp/dy) 0 ^^*2 


(27) 


5 

For not too small distances from the wall for not too small Pr 
numbers 1 « A^/n < Pr^A^/p, the integrand of equation (20) becomes ap- 
proximately equal to k. Since in the immediate neighborhood of the 
wall, or small values of q>, k * 0 and the integrand thereby vanishes, 
the latter can be replaced by k practically for the entire region of 

the friction layer (hence, the approximation formula t = 1 kdcp) . 

Jo 
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it is convenient to introduce a special heat-transfer parameter. A fur- 
ther question to be considered here, however, is whether the new dimen- 
sionless parameter is to be formed with u* or with U. 

The employment of U has the advantage that a good connectioh is 
made with older formulas . There will therefore be defined a "generalized 
heat- transfer coefficient" a* through the equation 


a* = 


10 

P0 c p 0 ® u 


(28) 


whereby the heat flow at the wall is referred to the maximum value of 
the temperature difference and of the velocity in the boundary layer. 


In place of equation (28) we can also write 


(28a) 


or 


a* 


Nu^ u cp 


m 


Pr 0 R e 


(28b) 


The expression Ru/(Pr Re) is frequently used in technical literature. 
Since and <p m lie in the neighborhood of 1 (figs. 7 ahd 9), a* is 

of the same order of magnitude as lfu/(Pr Re). 


Using equation (28) gives, fiom equation (26), 


a* 



(Pr* - l)d<p 
1 + Pr*A T /[i 



(29) 


This formula differs from equation (2l) , among other respects, in 
that here l/(l + Pr^A^/p) = q^/q is integrated over <p, whereas in 

equation (2l) the corresponding integration is made over y/r. Since 


f 


fr T a(y/r) 

and , on the other hand , — .. decreases with the distance from the 

<3 (y/r) 
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wall - for logarithmic velocity distribution 


is proportional to 


d<p 

a(y/r) 

l/y - it follows that the integrand of equation (29) decreases at a con- 
siderably greater rate with distance from the wall, than the correspond- 
ing integrand of formula (21) . 


For this reason, in using formula (29) only the relations in the 
immediate neighborhood of the wall need be known if the Prandtl number 
is not too small. It is, therefore, a particular advantage of formula 
(29) that the knowledge of the results of the universal wall friction 
law, in general, is sufficient for its- application. A. further advantage 
is that formula (31) for Pi* = 1 becomes particularly simple. 


P 1 S 

The magnitude e =s I kd<p occurring in formula (29) is only a 

JO 

small magnitude with regard to the main region of application of this 
formula (medium and high Prandtl numbers). Hence, in spite of the ex- 
isting gaps in our knowledge, the e-term can accurately be computed if 
the fl6w is sufficiently well known. For the case of the fully developed 
flows in a pipe and in a plan£ rectangular channel, e has been deter- 
mined by the author for various Pr and Re numbers (these values are 
given in fig. 8) . 

r 1 

Since the integrand k of the expression e = I kd<p differs ap- 

Jo 

preciably from zero only for the parts of the boundary layer that are at 
a large distance from the wall, the individuality of the friction layer 
under consideration is expressed precisely in the e-term. As may be 
seen from the smallness of the e-values, the individuality of the flow 
in turbulent heat transfer plays only a small part if the Pr and Re num- 
bers are not too small. 


6 


The term 



depends on the ratio (q/qo) / ( tt/tq) ( see e< l- 


(10)). This ratio is to a large extent determined by the boundary con- 
ditions at the wall (dq/dy)Q, or (dr/dy) q • In the case of pipe flow, 

(q/q^) / (t/tq) , and, therefore, k and s, is relatively large, since on 
account of the boundary conditions (di;/dy)Q < 0 and (dq/dy)o > 0 the 
distributions q/q^ and t/tq deviate relatively widely. In the case 
of the plane channel flow, q/qQ and t/tq no (Continued on next page.) 
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The principal problem of the turbulent heat transfer lies in the 

integral expression 

fore justifiable to introduce a special notation for this expression. 

It would be of advantage if the factor (Pr* - l) could first be taken 
outside the integral. Because of the local.. variability of Aq/A T and 

of the temperature dependent Pr, this unfortunately cannot, in general, 
be done. 


dtp of formula (29) . It is there- 
0 1 + Pr*Ar/M V . 


Accordingly, (Pr* - l) can be compensated only under the integral 
sign, and for this purpose it will be divided by the factor (Pr g - l) , 

where Prg denotes a mean value of the Prandtl number for the boundary 
layer near the wall. With regard to a later application of the universal 
wall friction law (see below) it is convenient, in addition, to repla ce 
the velocity <p by the universal velocity u/u#, where u* = V^o/PO 

and is the so-called shear stress velocity. The following magnitude b 
will therefore be defined as follows :' a 


j"iU/u* 


(Pr* - l)/(Pr g - 1) 
1 + Pr*Ar/(i 


d (u/u*) 


(30) 


When this expression is used, equation (3l) is obtained from equation (29) : 

\2 

(31) 


a* = 


(A q /A T ) m (u*/ur 

1 + (Pr g - l)b(u*/u) + 6 


longer differ so greatly, since In this case (dr/dy )q < 0 and 
(dq/dy)o =» 0. Hence, the e- values for the channel are smaller than the, 
corresponding values for the pipe, as may be seen, from figure 8. 


In the flow at the flat plate, the distributions of q/q^ and 
t /tq coincide approximately over the greater part of the <p -region, be- 
cause here we have (dv/dy)o = (dq/dy)o- Prom this it follows that the 
e-values for the case of the plate are still considerably smaller than 
the corresponding values of the plane channel flow. Hence, a general 
neglect of the e-term for the flow along the plate should be permissible. 


7a 

This magnitude, which in previous papers was denoted by a, has 
been redenoted by b in order to avoid confusion with the notation pre- 
scribed by the standard regulations (DIW 1341) for thermal diffusibility 
a. 


3871 
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In this representation of the turbulent heat transfer, the basic form of 
the old Prandtl equation (see eq. (41)) is expressed clearly. 


In order to obtain this customary form of the representation, the 
heat-transfer coefficient a.* would be formed with the maximum velocity 
U (see eq. (28)). There is still the possibility, however, of referring 
the heat-transfer coefficient to the shear-stress velocity u*, hence, 
of forming a dimensionless qQ/(poCpQ0u*) (as was previously indicated). 

We shall, therefore, consider briefly this heat-transfer coefficient. 


Since the identity relation 


^0 

PO^Po 011 ’ 


M, 

2 ^ 0 % 


holds, we. may write 


Nu ^u ^ U 
2Pr 0 n r *" a 


(32) 


(33) 


in place of equation (28) . With the aid of formulas (31) and (33) , the 
Susselt number can be directly determined without computing the tempera- 
ture -distribution. The practical application possibility of these for- 
mulas is, however, connected with the condition that the product Pr « 

does not exceed a certain limiting value (more details on this are given 
in the section on the problem of the heat-conducting layer. 


The validity of the formulas derived in the present section is de- 
pendent only on the single assumption that the heat flow and the momentum 
flow have approximately the same direction at right angles to the wall. 
This assumption is practically satisfied in all setups for which dis- 
continuous changes of the wall temperature are avoided. 

In order to apply the general formulas, specific assumptions must 
in addition be made on the magnitudes Ag/A T , A^/p, and t/tq, considered 

as "given". This will be done in the following section. 


In previous papers there is found in place of (u*/u) the so- 
called resistance coefficient £, which Is connected with u*/U by the 
following equation: (u*/u)^ = Tq/ ( pU^) =» 0.125 . 
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3. INTRODUCTION OF SPECIAL ASSUMPTIONS 
(a) Ratio of Exchange Magnitudes Aq/A^ 

Up to the year 1932 , in which Page and Falkner (ref. 5) observed in 
the wake of a heated cylinder the various possibilities of the propaga- 
tion of heat and momentum., the general idea has been held that the tur- 
bulence mechanism for the momentum transport and for the heat transport 
are identical. The tests carried out. by Fage and Falkner at the insti- 
gation of Taylor showed , however , that the heat is more strongly propa- 
gated than the momentum. This result supported the so-called vorticity 
transfer theory of Taylor (ref. 19), in which different mechanisms are 
postulated for the mass or heat transport and for the momentum 
propagation. 

The stronger propagating ability of heat as compared with the mo- 
mentum transport was also confirmed later in free jets. According to 
the investigations so far available, the ratio of the two exchange mag- 
nitudes in free turbulence is given by A^/A^ * 2. 

These results had as yet no 'effect on the theory of heat transfer 
in turbulent friction layers. In this connection, there were considered 
the velocity and temperature measurements of F. Elias in a turbulent 
airstream along a heated plate (ref. 3, 1930). From these tests a far 
reaching congruence was obtained of the temperature and velocity pro- 
files in friction layers . This result was looked upon as a confirma- 
tion of the previously held conception of the identity between momen- 
tum exchange and heat transport. 

This conclusion from the measurement results of Elias was , however , 
in error, as was shown by the author in the year 1940. From the agree- 
ment in the dimensionless <p and £ profiles (or from (d^/dcp)o >» l) , 

the agreement of A^ and A T holds only under the assumption that 

Pr =.l. In the tests of Elias, who used air as the flow medium, Pr was 
approximately 0.72. For this Prandtl number and for (d$/d<p)Q =» 1, it 

follows from equation (18) , however, that Aq/Ax « l/o.72 ® 1.40, that 
is, a value considerably greater than 1. 

Tests of Lorenz and Friedrichs (ref. 9) , who measured the distri- 
butions of the velocity and temperature of an air stream in a heated 
pipe, likewise led to a similar result. (Here the author obtained the 
value Aq/A T » 1.5 from the experimentally found value (d£/d<p)o * 0.97.) 

These tests show, therefore, that the heat exchange in friction 
layers is also greater than the momentum exchange. The ratio Aq/A^ * 2 
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observed in free turbulence is, however, not attained in friction layers. 
For this case we apparently have 

1 < (A q /A T ) m < Ag/Aj < 2 (34) 


This is supported by the following considerations : 

The fact that (Aq/A r ) m in friction layers is smaller than in free 

turbulence permits concluding that the lowering of this ratio is a con- 
sequence of the effect of the wall. Such an effect is naturally strong- 
er in the immediate neighborhood of the wall than at a point at a large 
distance from the wall. If, therefore, the ratio of the exchange magni- 
tudes is lowered through the effect of a wall, it follows that Aq/A T 

decreases with decreasing distance from the wall (or increases with in- 
creasing distance from the wall) . 

The value of the magnitude that A^/A t apparently approaches at a 

large distance from the wall is 2, which is observed in free turbulence. 
On the minimum value of Aq/A^, in the Immediate neighborhood of the 

wall, no experimental observations are as yet available for this region } 
therefore, a hypothesis must be assumed. The most probable assumption 
with regard to the neighborhood of the wall would be that only a single 
exchange mechanism holds at the wall, at which, therefore, (Aq/A^)^ = 1. 

This leads us back to the old hypothesis of the equality of the exchange 
mechanisms , but with the restriction that this identity of heat exchange 
and momentum exchange is as sinned only for the immediate neighborhood of 
the wall . 

There will now be considered the consequences that follow from this 
hypothesis. The magnitude Aq/Air appears in equations (31) and (30) , 
where it is contained in the magnitude Pr* = Aq/Ap Pr. Equation (30) 

will be considered first. Under the assumption of an explicit wall fric- 
tion law (for large Pr and Re numbers), the integrand of equation (30) 
exists only for the immediate neighborhood of the wall, for which there 
is to be set Aq * At* In equation (30) , therefore, Pr* can be replaced 
by Pr., This is of great use in the further considerations. 

The factor (Aq/A T ) m , appearing in the numerator of expression (31) , 

is represented by equation (23). This formula may be simplified by set- 
ting approximately cp = cpQ (on account of the weak temperature depend- 
ence of the specific heats , this simplification should always be permis- 
sible) . In place of equation (23) we can, therefore, write approximately 



d* 


At <3(y/r) 


y 

r 


(23a) 
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The factor (Aq/A T ) m , therefore, represents a mean value of the ratio 

of the exchange magnitudes over the entire region of- the friction layer. 
This mean value is to be formed directly over # or indirectly over y/r; 
the variable y/r is introduced because A^/a^ does not depend directly 

on but must be regarded as a function of y/r (or of r] r y/r ) . (For a 
preassigned Reynolds number, Tj r = ru*/v is a constant, as will be seen 
later . ) 

Since the variation of 7^ (y/r) is as yet unknown, only statements 

of a general character can be made at first with regard to the mean value 
of this ratio. For example, from equation (23) we may derive the result 
that the condition (34) is satisfied if A (J /a t has the limiting value 1 

at the wall and the value 2 at the edge of the friction layer. 

A further result is that the Prandtl number plays a part since ^ 
and d^/d (y/r) depend on this number. There may first be considered the 
relations for a large Pr where the temperature as is known is sufficient- 
ly well balanced over the flow cross section (fig. 4). The temperature 
gradient is, therefore, approximately zero over the cross section. In 
the layer near the wall, however, d0/d(y/r) reaches considerable values. 

On account of the small thickness of this layer, however, A^/A x may be 

treated as a constant , and according to the assumption may be set equal 
to 1. There follows, therefore, from equation (23a) 

(A q /A T ) m -► 1 for Pr -*■ - 

The value m thus attains the assumed minimum value at high 

Prandtl numbers. On the other hand, from equation (23a) it is seen that 
as this mean value deviates more strongly from 1 the sharper the temper- 
ature profile at relatively large distances from the wall where A q /A T 

approaches the assumed maximum value 2. According to our hypothesis, 
therefore, the smaller the Prandtl -number, the greater the values of 
(Aq/A^ja to be expected. 

An explanation of the questions raised here, also in particular the 
question of the effect of the Reynolds number on (A q /A T )m, can be pro- 
vided only by experiment. Of the test results thus far available, the 
measurements of the temperature and velocity distributions undoubtedly 
favor the previously presented view on the effect of the wall on the 
ratio of the exchange magnitudes. 


T/.fi9 
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Under the assumption that in the denominator of equation (18) we 
can set Pi* = Pr ■ and with account taken of a new determination of e 
(ref. 15) , there follows "both from the test results of Elias on the flat 
plate and from the results of Lorenz and Friedrichs (ref. 9) in the pipe 
according to equation (18) : 

(Aq/A T ) m « 1.3 for Pr =* 0.72 

According to the above considerations , this mean value , on account of 
the relatively small Prandtl number of the air, appears plausible as re- 
gards its order of magnitude. 


to 

A 


(b) Ratio of Friction Coefficients A^/p 

If we disregard the region of the boundary layer near the wall, the 
ratio of the friction coefficients A^p is an individual function of 
the distance from the wall for the case under consideration. 

For the fully developed flows through pipes and channels, A^/p may 
be represented approximately by the following function (ref. 16) : 

is- » * [o.s + ( Z /r) 2 ][x - (z/r)2] (35) 

where z/r = 1 - y/r is the distance from the axis of the pipe or chan- 
nel (A^/p must be symmetrical with respect to the axis), and x is the 

universal constant known from the mixing- length theory® which has an ap- 
proximate value of 0.4. The effect of Reynolds number is expressed 
through the factor Tj r = ru*y(> . The function (35) is plotted in fig- 
ure 1. 

As may be deduced from formula (35) or seen from figure 1, A T / (pT) r ) 

increases near the wall approximately linearly with the .distance from 
the wall with the slope 

A^/C^r) 

a X 

On account of q r = ru*/i , it follows that for A^/p 

A^/p = xt] (36) 



®See L. Prandtl: Fuhrer durch die Stromungslehre . Braunschweig, 

1949, pp. 118-119. 
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Near the wall, therefore, A^/p is determined hy the dimensionless 

distance from the wall t] = yu*/v, and the distance y/r no longer ap- 
pears. It should be borne in mind, however, that formula (36) holds 
only for the "fully turbulent region",' that is, where A T /[i is suffi- 
ciently large. The turbulent exchange vanishes, however, not only at 
T] - 0, but is actually already very small at a certain distance ti • 

With regard to these relations in the immediate neighborhood of the 
wall, the author has set up the following expression® (ref. 15) : 

At/m - *[t) - % ^g(T)A n )] (37) 

where t] n is a measure for the thickness of the viscous wall layer. 

The corresponding curve is shown in figure 2 for the values of the con- 
stants x = 0.4 and t} n « 11, derived from flow measurements. The dis- 
tance for which A r /n = 1 is T] - tji * 10.8 (thus approximately 
ni - %)• 


J 


Substituting expression (37) into the relation obtained from ex- 
pression (17) 


du/u* ^ jho 
dt} 1 + A^ii 


(38) 


there is obtained for t/tq^ the universal velocity distribution u/u*, 

which is likewise represented in figure 2. As may be seen from the po- 
sition of the plotted test points of the author (ref. 15), this computed 
velocity variation is experimentally confirmed at least for not too smal 1 
rj-values (for very small ^-distances the experiments are too inaccurate 
to serve as a test of the theory) . 


For extreme r|- values equation (37) goes over into the following 
function : 


A^/p = x(i} - ri n ) (for t) » T| n ) (37a) 

- ( x /3t|n)il 3 (for T) « T) n ) (37b) 

Equation (37a) differs from (36) only in that the surface of the wall ap- 
pears displaced toward the fluid by the thickness r| n of the viscous 
layer . 


^Where Tang X = f* ", 

e x + g-X 

■'•^In the neighborhood of the wall where only the Ti-distance (and not 
y/r) is in question, we can set t « tq- 
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Equation (10) should be corrected as follows: 



dx dy 
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Since from equation (36) or (37a) there follows the well-known loga- 
rithmic velocity profile (44), which was always observed for large t], 
the expression (37) corresponds to the experimental data for these values 

Of T) . 

With regard to the variation of A^Yp for small t), only supposi- 
tions can be made. The turbulent longitudinal oscillations u ! are to a 
first approximation proportional to the distance tj. From this the re- 
sult may be derived, for a two-dimensional oscillatory motion, with the 
aid of the continuity condition, t hat t he transverse oscillations v* are 
proportional to Since = pu T v ! , it follows from this considera- 
tion, if u T and v 1 * * * are correlated , that A ^ The expression 

(37) should, therefore, be physically questionable also with regard to 
its variation for small t) (eq. (37b)), 

In the application of equation (37), however, at least for very high 
Prandtl numbers, heat-transfer coefficients that were too high were com- 
puted (see more detailed discussion, sec. 5). From this it follows that 
the expression (37) for very small t|- distances gives too high values of 
A^/p. Probably the assumed correlation of the longitudinal and trans- 
verse fluctuations does not exist in the immediate neighborhood of the 
wall, and the power series for A^/p first begins with a term of the 
fifth degree. The variation of A^/p for the region T| < 6 should, 
therefore, be corrected and written 

(A c /n) x = 2.7 • 10~ 5 t) 5 (39) 

From r) > 9 on the original variation according to expression (37) is 
maintained (see table i). In figure 2 the corrected variation of A^p 
is shown by the dotted curve. 

This change of Arr/p has as a consequence a velocity profile (shown 
dotted) which hardly differs from the original profile (the deviation is 
at any rate considerably smaller than the scatter region of the test 
points) . Although the change made in A^/p for very small tj- values 
hardly has any hydrodynamic effect, the thermal effect of this correction 
at high Prandtl numbers is quite considerable , as will be seen from the 
considerations in section 5.* 


1 ^The third power of tj is the lowest power that can be considered. 

The expressions of K. Elser (ref. 4) and R. G. Deissler (ref. 2) which 

contradict it cannot, therefore, be corrected for the immediate neighbor- 
hood of the wall (here they give exchange values that are too high) . 

Remark made during the proof correction: 

In further investigations a new. expression (Continued on next page.) 
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(c) Computation of Temperature Distributions and Nusselt Numbers 


The formulas given above served for the computation of the tempera- 
ture profiles of turbulent flow in a pipe. This computation was made on 
the basis of equation (18) for different Prandtl numbers Pr (between 0 
»nr| 1000 ) , which were locally invariable and for the Reynolds number 
Re = 3 x 10 4 (or t) r = 800) . For all Prandtl numbers there was thus as- 
sumed a definite velocity profile <p , namely, the profile of the fully 
developed flow for which t/t q = 1 - y/r.^ The constants of the mate- 
rial Pr and Cp were assumed as constant . 


The magnitude A^/p occurring in equation (18) was determined for 
t) < 6 from equation (39), for ti < 30 from equation (37), and for 
r| > 50 from equation (35) . The temperature distributions of the first 
approximation (k = 0) determined by graphical integration served for the 
computation of the first approximations of q/qQ according to equation 
(14) . These distributions of the heat flux q/q q, similarly obtained 
by graphical integration, are plotted in figure 5.^ 3 


From q/qo there is obtained 1 + k = 


<iAlo 


and with the aid of 


1 - y/r’ 

this relation the d-profiles can he computed to a second approximation. 
These second approximations of which should not deviate much from 
the exact solutions, are plotted in figures 3 and 4. 


By the same numerical approximation procedure, there was also com- 
puted the temperature distribution of the laminar flow with the aid of 
equation (13) (the parabolic velocity profile was assumed here) . The 
profile obtained in this manner, shown dotted in figure 4, does not ap- 
preciably differ from that of L. Graetz (ref. 6) and W. Nusselt (ref. ll) 
computed by an exact method . 


\/n 53 *T| n [ TlAn ■ Tan e(TiAn) - V 3 Tang 3 (T}/r) n )] (39a) 

was found to give good results for x = 0.4 and ri n =* 7.15. The dis- 
tance from the wall for which - 1 is here found to be th ■ 10.7. 

For very small r]- values expression (39a) goes over asymptotically into 

A v /n = 0.2 XimOlAln) 5 - 3xlO- 5 T) 5 (39b) 

12 The numerical values of <p(y/r) for different Re numbers are 
given in. reference 15. 

-^The heat flux first increases with y/r because the cross section 
2rt (r-y) Ax, through which the radial heat flows, decreases more strongly 
with y/r than the radial heat flow itself. 
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With the temperature distributions there were simultaneously ob 


tained also the temperature gradients at the wall 


(«$•)<>' 


which accord- 


ing to formula (20) are equal to 0-5 Nui u . Also is obtained from 

the temperature computation (fig. 9). Table 2 gives the values of Nu^ 


r u> 


£ u , and Nu, and in figure 6 log Nu is plotted as a function of log Pr. 

As may be seen from this curve, the Nusselt number , with decreasing 
Prandtl number, asympotically approaches its limiting value for Pr = 0, 
which is indicated by a horizontal line (this value is about 5.22). The 
dotted line below it corresponds to the value Nu =* 3.77 of the laminar 
flow (according to the exact computation of Graetz and Nusselt, Nu « 3.66). 


Further, with the aid of equation (13) the Nusselt numbers were 
computed for vanishing Pr for various Re numbers (table 3) . The dif- 
ferences in the Nusselt numbers are conditioned only by the difference 
in the velocity profiles . ** 


4. THE PROBLEM OF HEAT-CONDUCTING LATER 

As regards the computation of b according to equation (30) , funda- 
mental difficulties arise through the fact that the Prandtl number (par- 
ticularly for viscous fluids) depends on the temperature which greatly 
varies precisely in the neighborhood of the wall. It is thus necessary 
to compute with a Prandtl number Pr* that depends on the distance t| 
even if it is assumed that A^ and A T are identical. 

Here we shall not, however, enter further into these complicated 
questions. The basic assumption of the complicated theory for variable 
"constants" of the material is the solution of equation (30) for un- 
changing values of these constants . In what follows , this case will now 
be considered. 


**Added remark during proof correction: 

At the conclusion of the above investigations a paper by R. N. Lyon 
appeared on the heat transfer of liquid metals (ref. 20) . Lyon computed 
the heat transfer on the basis of continuity considerations, the assump- 
tion, among others, being made that the temperature gradient in the flow 
direction dT/dx is independent of the distance from the wall. The de- 
rived formula is applied to the region of small Prandtl numbers , which 
are characteristic of liquid metals. The Nu values of Lyon T s theory 
are larger throughout than the corresponding values of the present study. 
For example, for the viscous flow, Lyon computes the value Nu = 4.36, 
whereas in the above computation the value 3.77 was obtained (and the 
exact value is 3.66). The heat-transfer coefficients obtained by Lyon 
are, on the average, smaller than his theoretical values. 
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Under the simplifying assumption that PTg = Pi* =■ const, equation 
(30) with account taken of equation (9) goes over into 


h 



du/u* 

1 + Pr A^p 



(30a) 


A rough approximation for b is obtained by starting from the earlier 
customary division of the friction layer into a "laminar" wall layer and 
a "fully turbulent" principal layer. In the laminar layer, which extends 
up to the distance rja, A-^/p = 0, whereas in the turbulent region 
(q > Tyb) we have A^/p -+■ • . In this idealization the integrand of equa- 
tion (30a) is 1 for the laminar layer and 0 for the turbulent region. 
There is therefore obtained from equation (30a) 

b = u a /u* =» <p a U/u* (41 ) 

where <p a = Ua/U is the dimensionless flow velocity at the bounding 
edge of the laminar layer. 

If this expression for b is substituted in equation (3l), A^/A T 
set equal to 1, and the e-term neglected, there is obtained 

a* = l " (^ ran< 3'tl equation) (4l) 


This, according to its physical content, is the old equation of L. 

Prandtl (ref. 12) for the turbulent heat transfer. 

A better approximation for b is obtained by dividing the friction 
layer into three zones : the viscous layer at the wall, the laminar- 

turbulent transition region, and the fully turbulent region, already 
discussed in the introduction (ref. 14). We shall not, however, repeat 
the discussion here, but in what follows we shall apply expressions (37) 
and (39), which represent A T /p as a continuous function of the distance 
T). In actuality, there exists no completely laminar layer, some turbu- 
lent exchange always exists even in the' immediate neighborhood of the 
wall. This exchange, hydrodynamically considered, may be very small; 
but this does not mean that the minimum exchange has only a slight ef- 
fect on the heat transfer. 

At high Prandtl numbers, in spite of "practically" viscous flow in 
the immediate neighborhood of the wall, there takes place a predominantly 
turbulent heat transfer, as may be deduced from equation (8). If, for 
example, we substitute A^/p =» 0.01 and Pr* = 1000 in this equation, 
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there Is obtained q^./q m =10. In this case, therefore, the turbulent 

heat transfer is ten times as large as the molecular heat transfer, al- 
though the turbulent exchange A T is negligibly small as compared with 

the viscosity p. From this example it is seen what errors are possible 
if the turbulent exchange in the immediate neighborhood of the wall is 
set exactly equal to zero . 

v 

This state of affairs makes it appear advisable to differentiate 
between a "predominantly viscous layer" and a "predominantly heat con- 
ducting layer". The viscosity is predominant in the region near the 
wall 0 < t) < r) 2 _ if tiie distance tj]_ is defined by the condition that 

and T m are equal there. On the other hand, the heat conduction is 
predominant in the region 0 < T) < ti 2 if the distance T]g is charac- 
terized by the fact that here q^ = q m , or Pr^A^/p = 1 (see eq. (8)). 
The defining equations for t)]_ and rjg &re therefore 

hi) = 1 (42) 

^ (rig) = l/Pr* (43) 

From these equations it follows that the layer thicknesses Tj-^ and 

rjg are equal only for Pr* = 1, and that 2 , and, therefore, rjg 

itself, decreases with increasing Prandtl number. At very high Prandtl 
numbers tj 2 li es deep within the predominantly viscous region near the 

wall. But with decreasing Prandtl number, the t| 2 -boundary is displaced 

far.'* into the predominantly turbulent region . 14 At extremely low Prandtl 
numbers, the heat conduction exceeds the turbulent heat transfer in the 
entire region of the friction layer. 

For the numerical determination of the magnitude b according to 
equation (30a), it is important to know at which point y^/r the limit 

rjg lies, so that it can then be decided which function for A T /p to 
substitute in equation (30a). Since y 2 /r = T] 2 /r) r , y 2 /r depends on ri r 
or on the Re number. 


14 If , in addition, the law for the wall friction holds in this re- 
gion, equation (43) may with the aid of equation (37a) be transformed into 


^2 



(43a) 


This equation holds, however, only for sufficiently high Reynolds numbers. 
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Universal values of -b, that is, not influenced by the Re number but 
determined only by the Prandtl number, are obtained if yg/r is so small 

(smaller than 0.05) that within the limits of the heat conducting wall 
layer A^-/ p or u/u* depends practically only on r| . For very small 
Prandtl numbers, for which T]g according to equation (43) is very large, 

this condition can be satisfied, however, only for very large r) r val- 
ues or extremely high Reynolds numbers . 


This universal b was now computed for different Prandtl numbers, 
with application of formula (37), which represents A^/p as a function 

only of t). This function serves to determine %/q =* l/(l + Pr Aj/p); 
cjjn/q was plotted, however, not against t) but against u/u* (fig. 10) . 


By graphical integration, b 



was then determined . 


These 


values of b are presented in table 4 and plotted in figure 11 as a 
function of log Pr. 


Before discussing the result of this computation, a more detailed 
consideration will be given to figure 10, in which the Q^/q curves are 

shown. In this figure is also drawn a dotted horizontal line 
q m /q =» l/2. Where this line intersects the q^/q curves of the individ- 
ual Prandtl lines, we therefore have c^j = q/2 or = q^. These 

points of intersection therefore characterize the respective limits of 
the region near the wall with predominant heat conduction. The abscissas 
of the intersection points represent the velocity Ug/u* at the bounding 

edge of the heat conducting layer for the individual Prandtl numbers. As 
may be seen, Ug/u* increases with decreasing Pr. 


For large values of ti for which the logarithmic velocity law 


u/u* 


i In n + const 
x 


(44) 


holds, c^/q may be computed in a simple fashion. From equations (30a), 
(36) , and (44) there follows : 

Qn/q = 0.5 + 0.5 Tangj^(^ - (45) 

This Qm/q distribution, holding for very high Re numbers and small 
Prandtl numbers,^ lies symmetrical to the point (Qn/q = 0.5; u/u* =* Ug/u*) , 

15j>or example, for Pr = 0.01 we have U2/U* = 19.4 (table 4). For 
this value there follows from equation (43a) or (Continued on next page.) 
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as may also be seen from figure 10. Hpnce, for large Re and small Fr 

»• 

-'U'/u* 

b « ^ du/u* « % (46) 

Jo q 

From this it follows further, if equation (44) is applied to u 2/ u *> 
equation (43a) (footnote 14) used for t| 2 , and the empirical values 
x = 0.4, const =s 5.5 and t| n = 11 are substituted: 

b » u 2 /u* = 5.75 log (11 + 2.5/Pr) +5.5 (Pr < O.l) (47) 

This formula shows that for very small Prandtl numbers, Ug/u* or b is 

approximately a linear function of -log Pr. This is also se£n from fig- 
ure 11. 

For medium and small Prandtl numbers, the relations are more com- 
plicated, The t^/q distribution no longer lies symmetrical to the 

point considered, and a certain difference arises between b and u 2 /u*. 

This difference is not large, however, (as is seen from table 4 and fig . 
ll) so that, for high Prandtl numbers also, we may set approximately 
b « u 2 /u*. 

Since for small u/u* values u/u* * t), we may also write for high 
Pr values : 


b « Ug/u* * 

For high Prandtl numbers, therefore, b is approximately equal to the 
thickness t \2 of the heat conducting layer. 


equation (44) (if we set x = 0.4) t) 2 261. This position of ti should 

lie very near the wall, that is, y 2 /r should be small since the univer- 
sal wall friction law here applied holds only under this assumption. If, 
in order to satisfy this condition, we set y 2 /r = 0.05, then there is 

obtained Ti r = 5220. To this value of T] r corresponds u/u* a 27.8 
and the Reynolds number Re = 2Ti r <p m u/u* = 2 .5x10^ • Furthermore, 
tj> 2 = u 2 /u = 0.70. For Pr = 0.001 the corresponding values are 
u 2 /u* = 25 (compare table 4), ti 2 = 2510, ri r = 50,200, u/u* =* 33.5, 

Re = 3*L0®, and q>2 = 0.75. 
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If we substitute the approximation formula b = ug/u* In equation 
(3l) , the denominator of this equation becomes 

1 + (Pr g - l)<p 2 + s 

where <pg = Ug/U denotes the dimensionless velocity at the outer edge*. 

of the heat conducting layer. By comparing this expression with the 
denominator of the Prandtl equation 1 + (Pr - l)<P a (eq. (4l)) where <p a 

is the velocity at the boundary of the laminar layer, the following re- 
sult is obtained: 

The new theory differs from the classical theory, among other re- 
spects, in the fact that in place of the velocity <p a at the boundary 

of the laminar' layer there enters a velocity bu^/U, which is approxi- 
mately equal to the velocity <Pg at the boundary of the predominantly 

heat conducting layer at the wall. 

Whereas the velocity <p a of the classical theory represents a 
purely hydrodynamic magnitude, the velocity. bu*/u » <Pg depends also on 
the Prandtl number, which is contained in the magnitude b «* ug/u*. The 

magnitude b is, however, determined exclusively through the Prandtl num- 
ber (that is, is independent of the Re number) , only if the Re number is 
so large that the heat conducting layer lies in the immediate neighbor- 
hood of the wall, where the universal friction law holds (i.e., all mag- 
nitudes depend only on t]) . This condition, which will be numerically 
formulated below, is to be observed particularly for small Prandtl 
numbers . 

Since at small Prandtl numbers and high Re numbers Pr * l/xqg 
(see footnote 15), we can under these conditions for the factor Pr 
appearing in equation (33) write 

Pr tk, « — — t— “ (48) 

xy 2 /r v ' 


If we now assume that the universal friction law still holds with suf- 
ficient accuracy in the region and in the neighborhood of yg < 0.05 r, 

there follows from equation (48) with x = 0.4: 

Pr r| r > 50 . . “ (49) 

If this condition is satisfied, the given. b- values may be substituted 
in equation (3l) . The factor bu^/u » <p 2 for such. small Pr (or such 

large b) then also remains always smaller than 1. 
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If in equation (49) T) r is further replaced by Re by means of the 
identity relation Re = 2T) r <p m U/u*, there is obtained as condition for 
Pr Re 


Pr Re > 100 %U/u* 

Using the numerical values in footnote 16, there is then obtained *■> 

Pr Re > 2500 (50) 

Accordingly, the universal b may be applied if, for example, for 
Pr ■ 0.1; Re > 2.5x10 , or for .Pr = 0.01; Re > 2.5x10 . 

5 . COMPARISON OF THEORY WITH EMPIRICAL RESULTS 

A good part of the extensive experimental material on* turbulent 
heat transfer has been evaluated by the Prandtl formula (4l) . It was 
immediately recognized that classical theory does not correctly repre- 
sent the experimental facts, in that the factor <P a ( which properly "should 

be a function only of the Re number) is also dependent on the Prandtl 
number . It was therefore considered as an essential object of experi- 
mental investigation to determine the dependence of the magnitude <p a 

on the parameters Pr and Re from tests . 

The results of investigations in this direction by various authors , 
in connection with turbulent pipe flow, have recently been evaluated by 
B. Koch (ref. 8). Almost all the authors make use of the original 
(Prandtl) formulation that <p a is proportional to Re - -'-*®. M. ten Bosch 

(ref. l) , however, assumes cp & = Re -0 * 1 . This formulation is undoubtedly 
the better one, because according to the very careful tests .of J. 
Nikuradse (ref. 10) u*/u can be represented for the usual region of 
Reynolds numbers by the equation 

u*/u = 0.125 Re -0 ' 1 (10 4 < Re < 10 6 ) (5l) 

From the series of empirical formulas, therefore, those of ten Bosch 
will be chosen and used for checking the theory. 

The formula of this author is [MCA note: See appendix.] 

(p„ = B pr~® *l®®Re - ® *-*- (ten Bosch, empirical) (52) 

where B is a coefficient which has the value 1.40 for the heating of 
the fluid and the value 1.12 for the cooling of the fluid. This relation 
is supported by numerous test results of various investigators > in par- 
ticular by the very careful measurements of G. Rohonczi (ref. 18). 
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If we set <p a = bu*/u in equatioff (52), there is obtained from 
equations (51) and (52) an empirical formula for b independent of Re: 

b a 8 B pr~0*185 (ten Bos^h, empirical) (53) 

Since ten Bosch developed equation (52) from measurements at medium and 
high Prandtl numbers, formula (53) derived from equation (52) must be 1 , 
compared with the theoretical b-values of the medium and large Prandtl 
numbers. For this purpose the b-values obtained through graphical inte- 
gration from table 4 will be represented by an approximation formula. 
This formula is 


b = 10 Pr -0 * 30 (1 < Pr < 200) (54) 

A comparison of equations (53) and (54) shows that for Prandtl num- 
bers in. the neighborhood of 1 the theoretical b satisfactorily agrees 
with the experimental b (the mean value of the factor B, which lies 
between 1.12 and 1.40, is 1.26). For high Prandtl numbers, however, the 
theoretical b-values are too small (for Pr = 100, for example, by the 
factor 0.6), that is, the theoretical heat- transfer coefficients are too 
high. 


This evident deficiency of the theory for high Prandtl numbers 
for which the heat conducting layer lies entirely within the viscous 
wall layer, is only to be explained by the circumstance that the ex- 
pression of equation (37) (which for the neighborhood of the wall goes 
over into eq. (37b)) , used for the computation of b, gives too high ex- 
change values for the predominantly viscous region. For this reason, 
as already stated in paragraph 3, section (b) , for the immediate neigh- 
borhood of the wall (t) < 6) there was set up expression (39), which 
gives considerably smaller values of A T /n (fig. 2). 

If these corrected values of A^/n are substituted in equation 

(30a) for the region t) < 6, there are obtained the corrected b-values, 
which are denoted by b k and which are likewise given in table 4 and 

figure 11. As is shown by figure 11, the dotted b k curve lies above 
the b-curve. At high Prandtl numbers b^ is considerably larger than 
b, while the difference at medium Prandtl numbers is insignificant (for 
Pr < 0.1 we have practically b k = b) . 

For the region of medium and high Prandtl numbers the previously 
given bjj. values may be represented by the following approximation 

formula : 


b k = 9.12 Pr-°*20 (l < Pr < 200) 


(55) 
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This formula agrees approximately with the empirical formula (53) . At 
any rate, the differences between equations (53) and (55) lie within the 
limits of the error range of the experimental results . 

A further correction of the assumed expression for the values of 
near the wall at the present time would also have no significance 

since the check of any theoretical fine details, in view of the uncer- 
tainty as to Ag/A^ and the uncertainties of the measured values, is 

quite impossible. In addition, it is to be borne in mind that the mag- 
nitude b (eq. (30)) is influenced by the temperature dependence of the 
Prandtl number (and, therefore, accurately speaking, is no universal con- 
stant)-'-® , and that this influence must be given theoretical treatment . 

The assumption by ten Bosch and other investigators that the factor 
B depends on the temperature gradient, in view of the structure of for- 
mula (30) , appears entirely plausible. Whether the problem, however, 
can be so simplified that two B-values are sufficient (for heating and 
cooling the fluid, respectively) is undecided. 


6. SUMMARY 

1. Under the assumption of a uniformly smooth wall, general formu- 

las are derived for the temperature distribution in turbulent friction 
layers and for the turbulent heat transfer. In these formulas the fol- 
lowing problematical magnitudes occur: the ratio of the turbulent to 

the viscous friction (A^/p) and the ratio of the exchange magnitudes for 

the heat and for the momentum (A^/a^) . 

2. In order to solve practical problems, concrete expressions are 
required for A^/A T and for A^/p. As regards the ratio of the exchange 

magnitudes, it was found from available tests that A„/A^ in boundary 

layers is lowered as compared with its value of about 2 in free turbu- 
lence . It is then assumed that Aq/A T with decreasing distance from 

the wall decreases from the value of about 2 at the outer edge of the 
boundary layer to the value 1 in the immediate neighborhood of the wall. 

3. For the ratio of the friction coefficients A^p in the various 

regions of the boundary layer, formulas are available which the author 
has developed in previous work on flow investigations (eqs . (35) and 
(37)). These formulas express the fact that the transition from the 
predominantly turbulent friction in the nucleus of the boundary layer to 


16 

The given h- values hold, 
ture differences . 


therefore, only for very small tempera- 
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the predominantly viscous friction in the immediate neighborhood of the 
wall takes place continuously. 

With the aid of these formulas for A T /|i, the temperature variation 
£ over the radius of a pipe (figs. 3 and 4;, the corresponding distri- 
bution of the heat flow q/q^ (fig* 5), and the Nusselt number (fig. 6) 

are computed by an already tested approximation method for a wide range 
of Prandtl numbers . 

4. The continuous transition from the predominantly turbulent to 
the predominantly molecular friction corresponds to a similar continu- 
ous transition from the predominantly turbulent to the predominantly 
molecular heat transfer in the neighborhood of the wall. The dimension- 
less thickness t] 2 of the wall layer with predominant heat conduction 

agrees, however, with the thickness of the predominantly viscous 

wall layer only for the Prandtl number Pr = 1. For high Prandtl num- 
bers the heat conducting layer lies deeply within the approximately vis- 
cous layer; whereas the heat conducting layer for very small Prandtl 
numbers may extend over the entire range of the boundary layer. 

5 . A computation of the heat transfer at the wall according to a 
general formula (eq. (3l)) is possible without preliminary computation 
of the temperature profile if the heat conducting layer is restricted to 
the region near the wall in which the universal wall friction law holds. 
This assumption is approximately satisfied if the product of the Prandtl 
number by the Reynolds number is greater than about 2500. 

6. In the general formula for turbulent heat transfer (eq. (31)), 
there occurs a factor bu*/U = <p 2 which corresponds to the factor <p & 

of the Prandtl formula (41) . Whereas the magnitude cp a of the classi- 
cal theory denotes the dimensionless velocity at the boundary of the 
fictitious "laminar layer" and thus represents a pure flow magnitude, 
cp 2 is approximately equal to the velocity at the boundary of the heat 

conduction layer; cp 2 is not a purely hydrodynamic quantity, since the 
factor b depends on the Prandtl number. ' 

For very large Prandtl numbers the magnitude b is approximately 
equal to the dimensionless thickness t) 2 of the predominantly heat con- 
ducting layer. 

7. The magnitude b which is independent of the Re number, de- 
pends on the heat- flow distribution q^q. This distribution (fig. 10) 

was computed using the expression (37) for A^/p. The b-values thus 
obtained (table 4 and fig. ll) agree for medium Pr numbers with the 
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experimentally obtained values as is shown "by the comparison with an 
empirical formula of ten Bosch. For high Prandtl numbers , however, too 
small b -values , that is, too large heat-transfer coefficients, are ob- 
tained by computation. From this it follows that expression (37) gives 
too high values of A T /^ for small distances from the wall (rj < 6) . 

For the region t| < 6 the values of Ap/p were therefore corrected 

in the direction of a still stronger decrease of the turbulent exchange 
at the wall (eq. (39)). The corrected values thus obtained bjj- (table 

4 and fig. ll) also for high Prandtl numbers, are quite well represented 
by the empirical .formula. 


A 




Aq 


B 

b 

C I 

d 

f 


k = qrro/ ( 10 T “ 1 
q 


T 

T - T 0 


SYMBOLS 

turbulent exchange magnitude, kg/m h 

for moment™ => A^ 

( dTN 

- c pAq ay ) 

factor used by M. ten Bosch 
factor used by H. Reichardt 

specific heat at constant pressure, kcal/kg (°C) 

diameter , m 

flow cross section, m^ 


heat flux (qQ at the wall; q m molecular heat flow), 
kcal/m^h 

radius of pipe or channel (flat sides) or thickness of 
the friction layer, m 

(absolute) temperature, °K 

temperature difference between flowing medium and wall 
( 0 maximum temperature difference) , °C 
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U 

u 


u m 

u* = Vw^o 

U, U' , V' 


X, y 


maximum value of the u- velocity, m/h 
flow velocity, m/h 

u-velocity averaged over flow cross section, m/h 
shear stress velocity, m/h 

flow velocities (u or u mean velocity; u* fluctua- 
ting velocity parallel to the wall; v’ fluctuating 
velocity at right angles to wall) , m/h 

coordinates (x parallel to wall; y distance from 
wall) , m 


z = r - y 


distance from pipe axis or from middle plane of 
channel, m 


Re " u ta d / v 


Reynolds number 


Pr =* n cp/X. 


Prandtl number 


Pr 


& 


Pr* = fa • Pr 
At 


Nu> a 


_*■ 
a ^ 


q 0 <3 

pc U0 

0 P0 


Prandtl numiber referred to a suitable mean temperature 
of boundary layer near wall 

generalized Prandtl number 
Nusselt number 

generalized heat-transfer coefficient 



coefficient 


£ resistance coefficient 

T) = yu*/v dimensionless distance from wall (t]^ thickness of 

predominantly viscous wall layer; thickness of 
predominantly heating wall layer; rj r a ru*/v) 
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® maximum temperature difference, °C 

4 s (T - T 0 }/® dimensionless temperature difference (£ m referred to 

mean temperature, to mean flow temperature) 


x constant of universal wall friction law 

X thermal conductivity, kcal/m h C 

p viscosity, kg/m h 

v == p/p kinematic viscosity, kg/m h 

p density of flow medium, kg/m? 

. du 

t shear stress (-to at wall; = p molecular shear 

stress; Tt turbulent shear stress), kp/m^ 

cp cp = u/U, <p m = Tim/U 


Subscripts : 
0 
1 
2 
a 

g 

k 

m 

m 

n 

r 

1 

t 

u 

T 


wall 

boundary of predominantly viscous layer 

boundary of predominantly heat conducting layer 

boundary of laminar layer 

boundary layer (near the wall) 

corrected magnitude 

molecular transfer 

mean value 

reference length 

pipe axis 

heat flow 

turbulent transfer 

flow medium 

momentum flow 
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APPENDIX 

DISCUSSION CONCERNING A CORRECTION TO THE TEXT 

The following extracts from an exchange of letters between Dr. 
Grigull and Dr. Reichardt are included in this translation at the request 
of Dr. Reichardt. ^ 

Dr. Grigull to Dr. Reichardt, March 4, 1952: 

I recently received your latest paper (Arch. ges. Warmetechn. ^ 
(1951) 129/42) which I read with great interest. I noted, however, that 
the empirical values of ten Bosch which were used for 'comparison refer 
to the ratio u a /u m (eq. (52)), whereas you refer the Prandtl equation 
(4l) to U. The values of. ten Bosch are, therefore, larger by approxi- 
mately the factor 1.2. If these values are reduced, the computed b- 
values will be found to lie between your original and the corrected bj^- 

values (table 4) . This means that your expression for the exchange 
values is more correct than appears from your comparison. 

I would appreciate it if you check my statements and let me know 
whether they are correct. I am particularly interested in this question, 
because I am preparing a new edition of my text book (Grober-Erk-Grigull, 
Grundgesetze der W&rmeiibertragung , 3rd revised edition. Springer -Ver lag 
1955) and would like to consider your theory in somewhat more detail 
since I consider it very suitable for arriving at a generally valid equa- 
tion for turbulent heat transfer. 

Dr. Reichardt to Dr. Grigull, March 13, 1952: 

You are quite right in your criticism. I have overlooked the fact 
that the <p of ten Bosch refers to the mean velocity instead of 

the maximum velocity U which I used . In my equation (53) there must , 
therefore, appear the £p m term. Thereby the deviation between the 
computed b-values and the experimental b-values at high Prandtl numbers 
becomes smaller. 


17 

Personal correspondence with T. M. Hallman, Lewis Flight 
Propulsion Laboratory, dated June 27, 1955. 
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TABLE 1 


r\ 

Ar/n 

(\/^ k 

0.2 

8. 8x10" 6 

9xlO“ 8 

• 4 

7 .IXIO -5 

3x10' 7 

. 6 

2.4X10 -4 

' 2.1X10 -6 

.8 

-4 

5.6x10 

-6 

8 . 8x10 

1 

.001 

-5 

2 . 7X10 

2 

.009 

8.6X10' 4 

3 

.030 

-3 

6.5X10 

4.5 

.093 

-2 

4.9X10 

6 

.209 

.193 

9 -* 

.630 

.630 

11 

1.05 

1.05 

12 

1.30 

1.30 

15 

2.14 

, 2.14 

18 

3.12 

3.12 

21 

4.19 

4.19 

24 

5.31 

5.31 

36 

10.0 

10.0 

48 

14.8 

14.8 

67 

22.8 

22.8 

100 

35.6 

35.6 
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TABLE 3 


Re 


*u 


3*io 4 

2.51 

0.496 

5.05 

3*10 4 

2.52 

.480 

5.24 

3X10 5 

2.55 

.470 

5.43 

3xio 6 

2.58 

.468 

5.50 


TABLE 4 


Pr 

0.001 

0.01 

0.1 

1 

10 

100 

1000 

^2 

2510 

261 

35.6 

10.8 

4.60 

2.18 

1.05 


2510 

261 

35.6 

10.8 

5.17 

3.24 

2.20 

u 2 /u* 

25.0 

19.4 

14.0 

8.84 

4.50 

2.07 

1.00 

(ug/u*)k 

25.0 

19.4 

14.0 

1 

8.81 

5.05 

3.25 

2.06 

■b 

25.1 

19.5 

14.1 

9.10 

1 

5.04 

2.49 

1.18 

*k 

25.1 

19.5 

14.1 

i 

9.19 

5.56 

3.47 

2.20 
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Figure 2. - Ratio of turbulence to viscous 
friction At/ \x and dimensionless velocity 
u/u* as function of the dimensionless 
distance from the wall, T| = 7\±*/v . 
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Figure 3. - Temperature distributions 
$ of turbulent flow in a pipe (sec- 
ond approximation) over the velocity 
<P for Re » 3X10 4 (n r * 800) for vari 
ous Prandtl numbers. 
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Figure 4. - Temperature distributions 
d of the turbulent flow in a pipe 
over the dimensionless distance from 
the wall y/r. Ep and Pr as in 
figure 3. 



Figure 5. - Distributions of dimension- 
less heat flux q/q 0 and shear stress 
t/t in a tube over the distance from 
the wa 11 y/r ( for Ee 3X10^ ) . 


























